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Abstract 



We compute thermal and quantum fluctuations in the background of a do- 
main wall in a scalar field theory at finite temperature using the exact scalar 
propagator in the subspace orthogonal to the wall's translational mode. The 
propagator makes it possible to calculate terms of any order in the semiclas- 
sical expansion of the partition function of the system. The leading term in 
the expansion corresponds to the fluctuation determinant, which we compute 
for arbitrary temperature in space dimensions 1, 2, and 3. Our results may 
be applied to the description of thermal scalar propagation in the presence of 
soliton defects (in polymers, magnetic materials, etc.) and interfaces which 
are characterized by kinklike profiles. They lead to predictions as to how 
classical free energy differences, surface tensions, and interface profiles are 
modified by fluctuations, allowing for comparison with both numerical and 
experimental data. They can also be used to estimate transition tempera- 
tures. Furthermore, the simple analytic form of the propagator may simplify 
existing calculations, and allow for more direct comparisons with data from 
scattering experiments. 
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I. INTRODUCTION 



There are many problems in physics which require quantization around classical back- 
grounds, such as Schwinger's QED calculation of vacuum polarization around constant elec- 
tromagnetic fields the semiclassical quantization of field theories 0, and the semiclassical 
analysis of the c/asszca/limit of quantum mechanics. In fact, the latter two characterize 
the semiclassical analysis as a two-way road: towards quantization, for the field theorists; 
towards the classical limit, for learning how classical physics descends from the quantum 
world. 

Taking the field theory lane of that two-way road, we have recently studied the problem 
of quantizing the simplest of field theories, one-dimensional quantum mechanics, whose 
(scalar) "field" lives in zero spatial dimension [^|-|^. Working in imaginary time, i.e., at 
finite temperature, and with smooth potentials, we were able to show that a full semiclassical 
series could be constructed from the knowledge of semiclassical propagators in the classical 
backgrounds of the solutions of the equations of motion, and gave a general prescription for 
deriving them from those backgrounds. Our construction was similar to the one proposed 
long ago by DeWitt-Morette |^ in real time, i.e., at zero temperature: it generalized to 
quantum statistical mechanics, where the number of classical solutions is drastically reduced, 
the results for ordinary quantum mechanics. 

This article applies some of our techniques to semiclassical quantization around a very 
specific background in a very specific scalar field theory, at finite temperature, in one, two, 
and three spatial dimensions. In that respect, it bears greater similarity to the Schwinger 
problem than to the semiclassical quantization of field theories. The latter would require 
finding all classical solutions that satisfy the finite temperature boundary conditions, and 
a general procedure for constructing semiclassical propagators in their backgrounds. As we 
have already stated, in quantum mechanics the general procedure does exist, and it is often 
possible to find all classical solutions In field theories, however, finding all classical 

solutions is already a difficult problem, and even if we find them all, no general procedure 
to obtain semiclassical propagators in their backgrounds is known. The technical reason for 
that is the presence of spatial dimensions, which transform the fiuctuation problem, given by 
an ordinary differential equation (ODE) in zero spatial dimension, into a partial differential 
equation (PDE) problem, whose solutions cannot be obtained from the mere knowledge of 
the background. Even finding the backgrounds involves, in general, a PDE, instead of the 
ODE of quantum mechanics. 

We have thus examined a specific potential of the yj'^-type, and confined our analysis to 
a (Euclidean) time-independent solution of the classical equations of motion depending on 
only one spatial coordinate, a domain wall whose classical profile and dynamics of fiuctua- 
tions could be exactly solved from ODE's. We were then able to construct a semiclassical 
propagator in the presence of that background, from which, in principle, we can generate any 
term of the semiclassical series around the domain wall. We note that even in one spatial 
dimension, where the domain wall has finite Euclidean action, and therefore contributes to 
the semiclassical approximation, we abdicate the idea of "solving" the theory semiclassically, 
because there are many other solutions whose semiclassical propagators cannot be obtained 
from a general recipe. In particular, solutions which only depend on Euclidean time do 
contribute: they are given by elliptic functions, similar to the ones in quantum mechanics 



2 



but whose fluctuation problems are much harder to solve. 
Our motivation for pursuing this problem, besides the technical appeal of obtaining 
thermal and quantum fluctuations at any temperature from an exact expression for the 
semiclassical propagator, was dictated by its potential physical applications. Indeed, for 
spatial dimension d = 1, our formalism can be adapted to describe fluctuations around 
defects that occur in one-dimensional systems, such as the solitons found in polymers (like 
polyacetylene) [|T^]; and for d = 2,3, to describe fluctuations around interfaces that separate 
two distinct phases |jlT|], a description that allows one to calculate how the surface tension, 
as well as the interface proflle are modifled by those fluctuations. It can also be used to 
approximately estimate transition temperatures |12[. Adopting this view, we are able to 
extract physically relevant results from the study oi d > 1, where the action of the domain 
wall diverges with the size of the system. 

Since the problem is treated at flnite temperature T, fluctuations incorporate thermal 
and quantum effects, which justifles the term thermal and quantum fluctuations employed 
thus far. The background we use, however, does not depend on T, as it is a "static" classical 
solution of the equations of motion, independent of the Euclidean time. That independence 
guarantees that it remains a classical solution at any temperature T, which is all we require 
to compute a semiclassical expansion order by order. As we incorporate thermal effects 
through fluctuations, our classical solution will no longer be a minimum of the (now T- 
dependent) Gibbs free energy (the effective action). Nevertheless, the new minimum may 
be obtained by extremizing the free energy with an ansatz that adds corrections to the 
original kink proflle. The procedure will be outlined in Section 0, and will proflt from the 
nice properties of that original proflle. Physically, in order for the procedure to be justifled, 
we must be dealing with situations where we start with a rather "massive" background, 
but whose fluctuations are "light", when both are compared to the temperature T. Thus, 
only the fluctuations would be influenced by the interaction with the heat bath, as in our 
calculations. This is very much what would be expected in the context of the so-called 
"adiabatic" approximations, where one can distinguish "heavier" and "lighter" degrees of 
freedom, as in the physics of polymers. Clearly, we have to assume that boundary conditons 
at spatial inflnity of the domain wall type are externally flxed, which forces us to go to the 
kink sector. 

We should stress that the problem of computing fluctuations around domain wall solu- 



tions has played an important role in the physics of interfaces in binary mixtures |TT[. In 
that case, however, temperature was introduced in a phenomenological way, through the 
"coarse-grained" Ginzburg-Landau free energy, whose mass parameter was assumed tem- 
perature dependent. In order to avoid confusion, we will call that temperature Tph. Our 
problem will reduce to that one if we also assume that the mass parameter of our Lagrangian 
depends on Tph, in a manner determined phenomenologically, and, furthermore, that our 
fluctuations are coupled to a heat bath at zero temperature, T = 0. Then, our methods can 



be used in the calculations that compute critical exponents to two-loop order ||T^, which 
obtain excellent agreement with both numerical and real experiments. Those calculations 
make use of the semiclassical propagator in the form of a series over eigenmodes. Our closed 
expression should simplify them, and should allow for an independent check. We, however, 
treat a different physical situation: one in which a microscopic fleld theory is forced to have 
two regions at distinct vacua. For that case, we compute how the interface separating those 
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regions is affected when a heat bath at T couples to its quantum fluctuations. In the polymer 
context, at c? = 1, that amounts to deriving how thermal phonons modify a soliton defect. 
For (i > 1, that can tell us how thermal neutral bosons affect an interface separating regions 
in different Bose-condensed phases. 

We should point out that a number of works have dealt with similar questions: the flnite 
temperature determinant has been computed long ago [|12],|14| , as a product over eigenvalues; 
Parnachev and Yaffe [|T3| have carefully studied the same problem at T = 0, paying special 

m 
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attention to the determinant calculation using a different method. Graham et al. 
investigated similar problems involving interfaces, resorting to the relation with phase shifts, 
rather than with the semiclassical propagator. In all cases, it seems possible to establish a 
precise relation to the approach we adopt. It should be noted that our approach includes 
a finite temperature T, and derives the exact semiclassical propagator which, together with 
a careful treament of collective coordinates, should allow for the computation of correlation 
functions beyond leading order. 

The structure of the article is as follows: in Section |I|, we briefly review how one de- 
rives the semiclassical expansion; Section [11 A| obtains an exact analytic expression for the 
semiclassical propagator, while Section |IIB| uses it to compute fluctuation determinants; in 
Section we describe the renormalization procedure, flrst by regularizing the calculation 
in Section |lll A , and then by removing its divergences in Section UIB ; in Section |V|, we dis- 
cuss the collective coordinates associated with the translational invariance of the solution; in 
Section 0, we show how surface tensions and interface proflles are modifled by fluctuations, 
and outline how the first higher order correction can be computed. We also mention how 
transition temperatures can be estimated from our results. We present our conclusions in 



Section VI 



II. THE SEMICLASSICAL EXPANSION 

The partition function for a self-interacting scalar field theory in contact with a thermal 
reservoir at temperature T can be written as {(3 = 1/T) 

Z{(3) = j[V^]e-^^^y\ (1) 
S[ip] = J^^ dr J d'^x C[ip] , (2) 

= lidr^r + li^^r + ^(^' - vir ■ (3) 

The integral / is to be performed over all y^'s that satisfy the boundary conditions (p{0, x) = 
(p{f3, x). The fields are rescaled so that the coupling constant A appears as an overall factor. 

We shall be interested in domain wall ansatze (p{x) which only depend on one spatial 
coordinate x = x^, which will be called longitudinal. As is well known [^, the equation of 
motion from @ is satisfied by a kinklike profile: 
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(p{x) = ±y3„ tanh 

where x denotes the position of the domain wall separating two spatial regions of dimension 
d. The remaining d — 1 transverse spatial coordinates of x define a vector r. The ± refers 
to kink and antikink, respectively. The classical action of the wall is proportional to the 
volume of the {d — l)-dimensional subspace orthogonal to the x-axis, which we will denote 
by Vd-i {Vo = 1): 

Sm = ^f3^lV,.,. (5) 

We expand the action around the domain wall configuration (p in order to derive the 
semiclassical series in that background: 

^{t, x) = ^{t, x) + A^/2^(r, x) , (6) 



(py{x - x) 
2VS 



(4) 



r/(0,x)=r7(/5,x) = 0. (7) 

As usual, the fluctuation r] has to vanish at r = and t = j3 because (p already satisfies the 
boundary condition. The expression for the partition function becomes 

Z{P) = e-'^^/' f[Dr]]e~^''^''\ (8) 

where both 5*2 and 5S depend functionally on (p and on t]: 



6S[^,V] = [drj d'x |^^r/3 + A^^l . (10) 

If we define = v^^/S, and use the explicit form of (p, the quadratic Lagrangian in can 
be written 

^2 = 1 [{drVf + (Vr/)2 + [M^ - \M\ed.\i)] r^^} , (11) 

with ^ = [ipvix — x)/2\/3] a dimensionless variable. 

We may generate a semiclassical series by expanding the e~^^ term, and using the 
quadratic propagator in the kink background to Wick-contract the various products of r] 
fields that will appear. Note that, when compared with perturbation theory, this expan- 
sion will involve different vertices, as well as a different propagator. The latter is already 
an infinite sum of one- loop perturbative diagrams [§,0. We now proceed to construct the 
propagator, a crucial step to obtain the series. 
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A. The quadratic propagator 



The quadratic propagator of our semiclassical expansion must satisfy 

j-a^- V2+[M2-^M2sech'(0]}G'(r,x;r',x') =5(r-r')5''(x-x') , (12) 

as well as the boundary conditions 

G(0, x; r', x') = G(/?, x; r', x ) = , (13) 

since the fluctuations vanish at r = 0, (3. 

The spectrum of the differential operator in (|T^ is that of a Posch- Teller potential, 
which is known exactly |T^. Many applications involving kinklike profiles |TT|- p!4| , p!8| have 
made use of that exphcit knowledge of the spectrum. Indeed, in some of them the quadratic 
propagator was expressed as a sum over the various eigenmodes of the operator, with the 
respective eigenvalues in the denominator. That, however, leads to very lengthy calculations 
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especially if one goes beyond one-loop [|T^. Therefore, it is useful to derive a dosed 
exact expression for the propagator which is actually the result of that sum. We can do it 
by realizing that, in transverse momentum space, the problem can be reduced to finding the 
propagator for an ordinary differential equation. 

Fourier transforming the transverse coordinates, as well as the Euclidean time, we obtain 
the ordinary differential equation 



i^-dl +ujI + \c' + [M^ - ^M^sech^O]] G{Kn, x, x') = 6{x - x') 



(14) 



where Un = 21112/ P, and we have defined the dimensionless vector k„ = -p-(co'„,k). As we 



(15) 



will see, G only depends on = |k„| = jj^^/iu^ + k^. The propagator G is given by 



G{r, x; r', x) = ^ / . w.^ G{Kn; x, x')e 

P n=-QO-' l^^J 



M V 



/^l^^aJ„(T-r')g^k■(r-r') 



In terms of the dimensionless propagator F = {M/2)G, and of the variable ^, defined in 
(|ll|), the equation reads 



(16) 



We can solve it by finding two linearly independent solutions, 0i(^) and 02(0; of the homo- 
geneous equation [|]-|^ 



{-9|+4 + «:^-6sech2(O}0.(O = O 
The solutions are obtained in Appendix 0: 



u 



l-u 



f)n/2 



1 - 



6 



-u + 
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hn + l (6„ + l)(6„ + 2) 



u 



(17) 



(18) 
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02 (m) 



l-U 



U 



fen/2 



12 



hn 



-u 



(6„-l)(6„-2) 



u 



(19) 



where m = (1 — tanh,^)/2 and 6„ = + k^. Expressions (plSf ) and (p!9|) are well-defined 
for Kn 7^ 0, i.e., 6„ > 2. We will return to this point later on. The propagator can now be 
constructed (see Appendix!^), yielding 



T{Kn\ u, u) = ^{0i(^) h{u) Q{u' -u) + (j)i{u') (j)2{u) Q{u - u)} 



(20) 



where G(x) is the Heaviside stepfunction. Note that the factors outside the brackets in ( p!8[ ) 
and (|T9|) will cancel in (20) whenever u = u' . 

The propagator we have just constructed can only be defined for k„ > 0, or 6„ > 2. 
As is well known P,p!7|, for = the fluctuation kernel in (|16]) admits a zero eigenmode 
associated to translational symmetry, i.e., a zero eigenvalue solution of (|17D that vanishes at 
^ ±cxD. Explicitly: 
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sech2[^(n)] = 2V?,u{l 



u 



(21) 



Consequently, we can only define a propagator in the subspace orthogonal to that spanned 
by (pl|). Thus, we subtract the contribution of that mode to obtain 



riQ{u)rio{u') 



(22) 



The divergent part of the limit as — > of ?(/€„;«,«') is exactly cancelled by the second 
term on the right-hand side of (p2|), so that !"(/€„;«,«') is well defined and orthogonal to 
the zero mode subspace in the limit ^ 0: 



lim r'(K„; u, u') = - u(l — u)u'(l — u') . 



(23) 



B. Determinants 

The leading term of the semiclassical expansion of the partition function around a domain 
wall background is obtained by setting 55* = in expression (H). The quadratic path-integral 
requires collective coordinates for the zero- mode subspace |TP| , |2D[ ]. As we will see in Section 
[rV| , the remaining integration, in the subspace orthogonal to that eigenmode, leads to a 
fluctuation determinant which will be computed from the quadratic semiclassical propagator 
for arbitrary temperature T and c? = 1, 2, 3. 

We again Fourier transform the fluctuation fields in all coordinates but the longitudinal 
one: 

Vir, x) = ^ E J Vi'^n, x) e^--^^^-^ . (24) 
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The quadratic part of the action may be written as 



S.-.Ejj^S.i'^.). (25) 



5'2(k„) = - / dxr]{Kn;x) K{Kn)r]{Kn;x) . (26) 

Z J—oo 

We denote by K{nn) the differential operator in (III), and formally define its determinant 



as 



Kl" = [detKM]-'/' =M j m e-^^(^") , (27) 

where the path integral sums over all fluctuations satisfying ±cxd) = 0, and A/" is a 

normalization constant. Denoting by \rjj) the eigenmodes (Posch- Teller modes, which do 
not depend on /«„), and by 7^^^ the eigenvalues of K{Kn), respectively, we have 

K-'M = t^-^, (28) 

= + 7j , (29) 

where j denotes both discrete and continuum Posch- Teller eigenstates, and jj = 7^0 are 
their respective eigenvalues. By convention, j = denotes the lowest of the eigenstates. As 
we have already remarked, 7gQ = 0, which implies Aq = 0. Therefore, when in the kq = 
subspace, we must use a primed determinant Aq which, by definition, is the product of all 
but the lowest eigenvalue. 

We may now relate ratios of primed determinants (i.e., with the lowest eigenvalue ex- 
cluded) to the semiclassical propagator of the previous section: 

ln^= r dsTTf'{V~s) , (30) 

where the trace is T'{y/s; ^,^). We are ultimately interested in ratios of determinants 

in the presence of the domain wall to free ones: for kq = 0, this is just given by the 
corresponding ratio for the Posch- Teller potential in one- dimensional quantum mechanics, 
ln(AQ/A^), which excludes the zero mode; for k„ 7^ 0, we may write 

A A' A-^ A' 

ln=?^ = ln«:^ + ln-^-ln-^ + ln^. (31) 

The first term in (^Tj) simply restores the lowest eigenvalue for k„ 7^ 0, whereas the second 
and the third can be calculated from the knowledge of propagators: the one in the domain 
wall background, and the free one. From the calculation outlined in Appendix O, we obtain 
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In- 



A' 



In 



^0 



1 / bn 

m 



bn + 1 



21n(6„ + 2) +ln48 . 



(32) 



We may follow to derive the result for kq = 0: 

In^ = - In 48 . 
^) into (|3TD , which is valid for k„ 7^ 0, we obtain 



(33) 



If we insert fB2f) and 



In 



A. 



Ar 



In 



+ In 



fen + 2 



(34) 



The above result includes the n = subspace of eigenmodes, but requires k„ 7^ 0, as the 
Kq = contribution has been explicitly summed to cancel the In 48 term of (^2|). 

Formally, the logarithm of the complete (primed) determinant is the sum over all possible 

K„'s > 0) of my. 



A' f df^^^k °° Ak 

In -TTT = Vd-i / -, — 73— r In 



(35) 



Expression (|35D , as it stands, is meaningless: the sum over n diverges logarithmically, 
and, for > 1, the integral over k introduces extra divergences as well. One must undertake 
a renormalization process in order to arrive at well-defined results. In the next section, we 
shall regularize, and then renormalize the calculation to obtain a final answer for any T and 
for d = 1, 2, 3, where the theory is renormalizable. 



III. THE RENORMALIZATION PROCEDURE 
A. Regularization 



We shall concentrate our attention on (pSl). We begin by examining the behavior of the 
sum S over Matsubara frequencies defined as 

N ^ 

E l^A^^' (36) 

n=-N Kn 



where iV is a cutoff. If we introduce a dimensionless function = (/3/27r)-\/M2 + k"^, the 
sum becomes 

^ + In ^ . . (37) 

/n2 + c| + |co ^n^ + cl + c^l 

S may be rewritten in a convenient way by resorting to expression ( |C5| ) of Appendix y. 



obtained from Plana's formula |21|, which is commonly used in the description of the Casimir 



effect [22|. For Kn 7^ 0, it is given as the difference of two terms: S = Sa — Sb, where 
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'5a = EE2e*'^- 

j = l (7=0 



N 



dx In 



+ cl 



(38) 



e^^y - 1 



vr — arctan 



JCo 



(39) 



The integral / in the brackets of (|38|) can be done exactly [see (|D3| ) in Appendix 0. For 
iV large, if we neglect terms of order replace with /5A/27r, and restore dimensional 
quantities, expression (|37D can be split into a zero temperature part Sq, 



5„ = ^a.c=inf££UM.,. + V„e,i4iL 

TT VCfc/ TT V 4 V2Cfc 



and a temperature dependent part St 
Sndy 



TT 



-(3M 



(40) 



arctan 



Co 



+ arctan 



2Vy^ - 4 
Co 



+ 4 In ( 1 - e 



■ (41) 



In both expressions, the first two terms represent the contribution of the (Posch- Teller) 
bound states in the longitudinal direction, while the last ones correspond to the continuum. 
For (i = 1, /c = 0, so that the first term of fEOl) vanishes. 



B. Renormalization 

The continuum contribution (last term) to the zero temperature part Sq contains the 
expected logarithmic divergence. Indeed, at zero temperature we should replace the sum 
(P^) by an integral: 



S 



N 

E In 

n=-N 



— In -r-F 

271 Af 



(42) 



where is a continuous variable that replaces and k, = -g(c(j,k). Furthermore, we may 
write 



In^ = Trln [l - («:)[/ 



(43) 



where 1 is the identity operator, tp^K,;^,^') is given by (|BTD, with 6„ replaced by 6 = 
V4 + K^, f/(0 = ~6 sech^(^) is the interaction term in ([T6|), and the trace is defined in (pOD. 
If we expand the logarithm in (|43|) , we shall have a series of one-loop terms: 



Aft °° 1 
In ^r-F = — — Tr 
tin 



(44) 
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Inserting (|4^) into (|42|), the n = 1 term Ti can be readily computed (change the integration 
variable from ^ to u) and identified with the divergent term of (^Oj): 



Ii = -P rd^Tp{K;^,0 f/(0 = - -/3Mln . (45) 

In order to renormalize the theory, we add a counterterm to the Lagrangian which will absorb 
that divergent term, and impose renormalization conditions to fix the finite remainder. For 
simplicity, we will adopt renormalization conditions at zero momentum k. In the present 
case, it suffices to add a mass-type counterterm: 

m^m-lc,{v'-vl), (46) 

and to require that the one-point one-particle irreducible (IPI) Green function at zero mo- 
menta be given by: 

f«(0)=0. (47) 

In order to satisfy (|47|) , we must subtract the tadpole graph in the kink background, which 
is also given by (|^) because the free propagator at coincident points can be factored out of 
the ^ integration. We then obtain the renormalized expression 



S^ = ^ arcsin (^) + ^Jk^ + ^M^ arcsin ( ^) - -/3M + St • (48) 

For d = 1, k = 0, Ck = Co, and there are no additional integrals. In that case, we may verify 
that the temperature independent part of (|i8| ) coincides with the result quoted in equation 



(3.29) of the first reference in [|1^. Indeed, from our calculation we may easily obtain the 



difference between the energy of the kink and that of the vacuum, i.e., the kink mass: 

1 1 
(3\ A 2 



M,,„, = ^ f ^ + ^^kI , (49) 



which at zero temperature yields: 

2M3 M f 1 3 \ , . 

The replacements M = \/2mraj, and A = GKaj, needed to account for different definitions 



of mass and coupling, show that our results coincide with Rajaraman's [1^ for d = 1. We 
note that ( P^j) and ( PB| ) had already made use of that coincidence for the case d = 0. 

For d = 2, and d = 3, we still have to integrate over k. This last integral will generate 
additional divergences in c? = 3 dimensions, which are contained in the n = 2 term of 
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1 



The calculation of Appendix |^ leads to 
J2 = -36 ' ^ 



+ 



1 



262(6 + 1)2 6(6+1) 



+ 0[l/(6+l) 



(52) 



(53) 



Since b = v4 + /?, the second term in (|53|) dominates in the ultraviolet. 

We now analyze expression (^Tj) for each dimension: i) in c? = 1, the integral over u goes 
like / uj~^ duo ~ A^^^ which vanishes as A — 00; ii) in c? = 2, the integrals over uj and k lead 
to I2 ~ A^^, which also vanishes in that limit; iii) in c? = 3, however, the integrals over uj and 
k lead to I2 ~ In A, which diverges in that limit. Thus, in c? = 3 an additional subtraction 
is required. The leading ultraviolet behavior of (^Tj) comes from the second term in 



X 



(as) 



d^k 



/9 



°° duj 



-6{M/2f 



-00 27T 



yAJ^+W+W + (M/2) 



(54) 



We integrate the momenta k over a sphere of radius A to obtain |23 

3 



(as) 



167r' 



V2(3M'' ln(A7M2) + 0(1) 



(55) 



This is to be compared with the result of integrating the finite terms of ( ^OD over momenta 
k (the divergent term was already cancelled by the firts counterterm) . Integration by parts 
leads to elementary integrals; expanding the result for large A, and neglecting terms of order 
l/A^, we obtain 



Vo 



_fk_ 



So 



V2pM^ ln(AVM^ 



V3 
487r 



V2PM^ + 0(1/A2) . 



(56) 



As expected, the leading ultraviolet behavior is the same, which shows that the n = 2 
term of ( ^^ is indeed responsible for the new divergence. The renormalization procedure 
requires that we include an additional counterterm in the Lagrangian, and use an additional 
renormalization condition to fix finite parts in d = 3 (see Appendix |F]). The Lagrangian 
becomes: 



\c2i^' 



'Vv) 



Ci and 6*2 are counterterms that we fix by imposing (0), as well as: 



fg)(0,0) 



(57) 



(58) 



The contribution of the additional counterterm to the effective action in the kink sector is 
given by: 



d^k r°° dw 



(27r)= 



I 



27r (w2 + F + M2)2 



4 Jo 



dr 



2\2 



(59) 
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where the first bracket is identified as the one-loop graph with two external legs at zero 
momenta (integrated over w in the interval [—00, 00], and over k in the two-sphere of radius 
A, as before), while the second is easily related to the classical kink action. A simple 
calculation yields: 



16tt' 



V2(3M^ ln(AVM2) + 



(60) 



If we subtract the result of ( |60|) from (|56|) , and take the limit A — >■ 00, we finally arrive at 



In 



167r2 



1 + 



TT 



3^3 



(Pk 



St ■ 



(61) 



We should point out that expressions for the finite temperature determinant were previ- 
ously obtained from the explicit knowlegde of the Posch- Teller eigenvalues, either by sum- 



ming over Matsubara frequencies, and then over those eigenvalues [0 , or by integrating over 
transverse momenta in arbitrary dimension first, making use of dimensional regularization, 



and then summing over Matsubara frequencies [O . 



IV. COLLECTIVE COORDINATES 

The domain wall solution breaks translational symmetry along the longitudinal x-axis, 
as it depends on a parameter x that characterizes its position. In order to restore that 
symmetry of the theory, we resort to the method of collective coordinates to sum over 



all possible values of x. In = 1, there is nothing else to restore. In ci > 1, where 
the wall also breaks rotational invariance by singling out a normal direction, consistently 
with our previous assumptions, the longitudinal direction is externally fixed by the forces 
that implement the wall boundary conditions at spatial infinity. Therefore, here as well, 
longitudinal translations are the only degrees of freedom to consider. 

The existence of longitudinal translation zero eigenmodes of the fiuctuation operator is 
a clear signal that the symmetry has to be restored. Physically, it means that it costs no 
energy to go from one solution to a longitudinally translated one. Thus, in the usual manner, 
in a region of longitudinal length L, we introduce in the partition function the identity 

df'x TjQ^x — x) [^{t, X, r) — 0{x — x)]j J[ip] = 1 , (62) 

where the Jacobian is given by 

jU=X-y^ dr /rf'^x(^(r,x,r)^r/o(x-x) = dr [ d'^x 7]o{x - x)^ . (63) 

JO J ox Jo J ox 

Equation ( |62D imposes orthogonality to the normalized zero mode r]o, and integrates over 
longitudinal positions. The second equality in ( p3| ) comes from an integration by parts. We 
note that the zero mode is already orthogonal to the domain wall profile since r^o oc d(p/dx, 
which implies 
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dxTjQ^x — x)(p{x — x) oc 



0. 



(64) 



Therefore, we may omit the (p term in (0). We insert the identity in the expression for the 
partition function (|l|), and change the order of integration to obtain 



(65) 



We perform the semiclassical expansion around the domain wall according to (^. Since the 
action can be written as 







dr I d'^xl^ 
\ ox 



(66) 



the normalized zero eigenmode is simply rjo = S[0] {d(p/dx). Using that, and its orthog- 
onality to (p, the expansion yields 



Z{(3) 



J ^e-^[^/^ j[Vi^e-'Hl^^f^dT J d'xvov) A^.v] E^^l"^^)"' (67) 



with the functional integral being performed over fluctuations that vanish at r = 0, /3, and 
the expanded Jacobian given by 



'S\g 



1/2 



dr I d'x 

ox 



'S\g 



1/2 



S[0Y^'^ dr j d'^xV'['p]ri . (68) 



The leading term in (pTf) corresponds to the quadratic approximation Z2 to the partition 
function in the domain wall background, whose ratio to the free one is given by 



Z2W) _ 


1 ^ g-S[^]/A 




-1/2 




Zp J 






R 





1/2 



(69) 



It involves the primed determinant calculated previously, since the delta function restricts 
us to the subspace orthogonal to the zero eigenmode, has a factor {S[0\/Xy^'^ coming from 
the Jacobian, and a factor of (27r) from the functional measure which is not cancelled 
by the corresponding value for the free determinant in the zero mode subspace. 



V. PHYSICAL APPLICATIONS AND HIGHER ORDERS 

The result obtained in the preceding section has an immediate physical application in 
the calculation of the interface tension in d = 2 and d = 3. If our system is confined to a 
longitudinal length L, and a transverse cross-section Vd-i-, the free energy difference per unit 
cross-section between the situation with and without the domain wall defines the interface 
(surface) tension |]TT]| 



cr = — lim — — — lim In ( — — ) , (70) 



Vd-i^oo j3Vd-l L^oo \Z 
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where Z denotes the partition function in the presence of the wall, and Zp the free partition 
function, without the domain wall. The lowest order result (|69D , with the integral over x 
taken from —L/2 to L/2 gives 



a 



lim 



1 f S[0\ 




(71) 



where we have left out terms vanishing as InVd-i/Vd-i or faster. Another way of deriving 
the interface tension is through its relationship with the energy splitting between ground 
and first excited state of the field theory, when its vacuum degeneracy is broken by the finite 



volume [£4|. The tunneling between the previously degenerate states can be computed using 
domain walls p5| , |26[] . This calculation was performed |^ using a dilute gas of kinks and 
antikinks [^, and yields 



5[?]/A 


"A'' 


-1/2 






_Ai. 


R 


V27rA ; 



1/2 



(72) 



an expression that comes from the exponentiation of (pOD , which results from summing over 
kinks and antikinks. Our results incorporate finite temperature effects into the expressions 



obtained in Ref. ||2^. Again, we stress that our starting point is a field theory that has no 
phenomenological temperature dependence in its Lagrangian, unlike treatments that start 
from a phenomenological Ginzburg-Landau coarse grained free energy ||lT|,|l3|. Clearly, for 
d = 1, a case that can be applied to the study of soliton backgrounds in one-dimensional 
polymers, we would be comparing the free energy difference with and without the soliton 
background (Vq = 1), and finding how it is affected by thermal phonons. 

The result obtained for the surface tension also provides us with an indirect determination 
of the temperature of the second order transition that will restore the —>■ —ip symmentry. 
Strictly speaking, it determines the so-called "percolation temperature", the temperature 
at which the surface tension vanishes, which coincides with the critical temperature in d = 
3, although it is smaller in d = 2 In fact, the statistical mechanics literature just 

compares free energies with and without domain wall boundary conditions, and establishes 
that the signal for the symmetry restoring transition is given when these two are equal. That 
means that the system will prefer to form walls separating regions at different vacua, whose 
condensation leads to a vanishing order paramenter. Our calculation is just the semiclassical 



version of that. That estimate has already been made long ago [0, using a different method, 
as already explained. 

Another physical information we can extract from our calculations is how the soliton 
{d = 1), or interface {d = 2,3) profiles are modified by thermal and quantum fluctuations. 
In lowest order, we may use the one-loop expression for the effective action. 



Am = ^ + lln{A[{cp)]) 



(73) 



where {<f{^)) is the expectation value of the field, and A is the Legendre transform of the free 
energy. We should look for extrema of (^), so we functionally differentiate with respect to 
((/^(x)). The resulting equation will just be the classical equation of motion plus a one-loop 
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correction which involves the propagator in the ((/^(x)) background. Neglecting the one- loop 
part, we know that the domain wall profile is a solution. Thus, the ansatz 

(^(x)) = ^(x) + X'/'5{^i^)) , (74) 

leads to an equation for the correction term that involves the propagator in the domain 
wall background, for which we have an explicit expression. As a result, we will see how the 
thermal and quantum fluctuations change the profile, and should be able to compare that 
prediction with scattering data from radiation (of appropriate wavelength) scattered off the 
domain wall. 

We close this section with comments about higher order corrections. In our construction, 
these are generated by expanding the e~^^ term in (^, and using the quadratic propagator 
in the domain wall background to Wick-contract the various products of i] fields that will 
appear. There is also an additional vertex that comes from the second term of the Jacobian 
in (pB|), which is linear in 7], and proportional to A^^^. The first correction to the lowest order 
result for the partition functions involves graphs which were computed in the second article 
of Ref. [r^l, and represent a two- loop correction. Those same graphs can now be calculated 



with our propagator either at T = or at finite temparature (r integrals going from to /5). 
In order to check and extend those results, one has to use the same renormalization scheme, 
conveniently chosen to facilitate comparisons with Monte Carlo data [|13|. 

Numerical calculations will be required in applying even our lowest order results, as 
well as to allow for comparisons with experimental or Monte Carlo data. The two-loop 
calculation just mentioned will also require a numerical effort. We plan to undertake such 
efforts in the near future. 



VI. CONCLUSIONS 

We have shown how one can systematically compute quantum fluctuations around a 
domain wall background which is temperature independent, but whose fluctuations interact 
with a heat bath at temperature T. The fluctuation determinant, the leading term in our 
semi classical expansion, already makes the effective action T-dependent, and one can obtain 
an extremum for it that corresponds to a dressed T-dependent version of the original domain 
wall. As higher orders are included, that dependence will obviously change. The determinant 
also leads to the calculation of the surface tension, which again will depend on T through the 



fluctuations. In both cases, one can follow the computations reviewed in Ref. , omitting 
their Tph- dependence in the paramenters of the Lagrangian, and replacing the determinants 
with our T-dependent ones. Also, the temperature T at which the surface tension vanishes 
can be used to estimate the phase transition to the unbroken phase [|T2 |. 



As we have already mentioned, results for = 1 can be applied to describe how thermal 
phonons affect soliton backgrounds in polymer physics (examples involving magnetic mate- 
rials are also potential applications), while those for d = 2,3 are applicable to Bose systems 
which are forced to condense in different phases in two different regions separated by a 
domain wall. The simple analytic form of the propagator in the domain wall background 
is quite useful in calculations, and should lead to the determination of physical quantities 
measured in scattering experiments. 
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In principle, the methods we use could be applied to other theories, as long as the back- 
ground solutions are sufficiently symmetric so as to depend on only one spatial variable (a 
radial one, for instance). The basic ingredient required is the knowledge, exact or approxi- 
mate, of two linearly independent solutions of the fluctuation ODE at zero eigenvalue, from 
which we construct the semiclassical propagator. Theories with solitonic (kinks, vortices, 
monopoles, etc.) or instantonic backgrounds often have such a property. 

We hope, in the very near future, to present detailed numerical calculations of the various 
quantities mentioned in Section 0. From them, we expect to establish how realistic our 
description is, by confronting it with available experimental data. 
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APPENDIX A: 

In this Appendix, we compute the solutions to (pT]), and use them to construct the 
propagator (pO|). The standard substitution 0(^) = e~"^ cosh~''"(^)F(^), with a = 



and 6„ = y4 + k^, and the change of variable u = (1 — tanh^)/2 lead to a hypergeometric 
equation for F: 

d'^F dF 
m(1 + [{bn + 1) - 2(6„ + l)u]— + [6 - 6„(6„ + 1)]F = . (Al) 

The general solution is the linear combination 

F{u) = ci 2Fi{bn - 2, 6„ + 3; 1 + 6„; u) + m"'" 2i^i(3, -2; 1 - 6„; u) , (A2) 
where 2Fi{A,B] C;u) is the hypergeometric function. Using the identity 

2Fi(A, B; C; u) = (1 - m)^-^-^Fi(C - A,C - B;C;u) , (A3) 

we obtain 

F{u) = c^{l-u)-'"2Fi{3,-2;l + b^;u) + C2U''"2Fii3,-2;l-bn;u) . (A4) 



Both series terminate, and we are led to the solutions in (|T8D and (|T9]). The propagator can 
now be constructed from the functions 



W{0 = 01 (002(0 - 0'i (002(0 • (A6) 
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It is given by [^-|^ 



where ^<(>) = min(max)(,^, ,^'). In terms of the variable m, noting that ^ > ^' ^ u < u', 
^ ^ oo ^ u ^ 0, and ^ —>■ —oo =^ m — >• 1, we obtain 

rK;M,M) = . (A8) 

The propagator does vanish when u,u' = 0, 1. Also, W{u) is a constant, easy to compute 
at M = 0; we find W{u) = 26„. We can then write expression (^). 



APPENDIX B: 



In this Appendix, we evaluate ( p2D from the propagators. The free one, which can also 
be obtained from the methods described previously, is given by 



-6n|5-ri 



(Bl) 



Changing variables to 6 = a/4 + s and u = {1 — tanh^)/2, the second term in (pl|) becomes 
A' 



In 



A'n 



2 J0U(1 — U) 



U,U] 



(B2) 



The expression inside the brackets is easily computed from the propagators, yielding 

6m(1-m) 12(^2 + 26 + 3)m2(i 

6(62 _ 1) 6(6 + 2)(62-l) ■ 

Doing the integrals in (B^), we arrive at 



(B3) 



APPENDIX C: 



In this Appendix, we shall make use of Plana's formula [pT] , pO| to compute the sum 
which can be decomposed into sums of the form 



N 



n=-N 



5(r,e^""s)= ln(Vn2 + r2 + e^""s) , (CI) 



with r > s > 0, and a = 0, 1. 

Plana's formula can be derived from Cauchy's theorem applied to two contours in the 
complex plane |21|]: one that runs counterclockwise in the upper half-plane, from {x,y) = 
(ni,ni + ioo) to the real axis along a perpendicular, then along the real axis, and finally 
perpendicularly away from it to (n2,n2 + ioo), avoiding the integers on the real segment 
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(ni,n2) (semicircling them), as well as the comers {ni,0) and (^2,0) (with 7r/2 arcs); the 
other is its mirror reflection on the real axis. For a function (j){z) analytic and bounded for 
rii < Re (2;) < n2, we have 



^ n=ni+l ^ 



n2 1 roo 

dx (j)(x) H — / dy 
I Jo 



ni 



Q2ny _ I 



(C2) 



where ^{n,y) = (f){n + iy) — 0(n — iy). Deflning (f)[z) = In (^V z"^ + + e^'^'^s^, we may 
rewrite (ICll) as 



S{r, e^""s) = 2 ^ In (Vn2 + + e^""s) - In (r + e*""s 



n=0 



(C3) 



and use (|C2|) with rii = 0+ (in order to avoid the square-root cut at Re{z) = 0) and n2 = N. 
The contribution to the second integral of ( |C^ ) involving ^{N, y) can be expanded for large 
N , and shown to behave as The one involving $(0^, y) has to be split into two pieces: 
\y\ > r and \y\ < r. The latter piece yields zero, whereas the former uses 



(j){0+±zy) = -\n[y'-{r^-s')]±z 



an + e*"" arctan f 



(C4) 



Neglecting terms that vanish as 00, we flnally arrive at 



S{r, e^""s) = 2J^ dx In [Vx^ + + e*""s) + In [VN^ + + e^""s 



-4 



dy 



r e 



(XTT + e^'^^ arctan 



(C5) 



where r > s > 0, and a = 0, 1. It is easy to verify that ( |C5| ) correctly reproduces the results 
for S{r,0) and for the sum Ei=o'5(r,e*'"^s), which can be obtained straightforwardly. 



APPENDIX D: 

In this Appendix, we calculate the integral / that appears inside the brackets of 

rN 

I{r, e'^^s) = dx In (Vx^ + + e*""s) . 

Changing variable to q{x) = \/ x"^ + + e'-'^'^s, we integrate by parts to obtain 

fiW dq r : : 

Ja O) a ^ 



q{N) 
9(0) 



{q — e'^'^'^s) — In g 
The last integral can be found in [^. Finally: 

J = AT In f VN^ + r2 + e'^'^s) - N + e'"^ s 



In 



< arcsm 



.2 ^ gia7r^^jY2 ^ ^2 



VA^2 _^ J.2 gi^TTg 



TT 

2 



(Dl) 



(D2) 



(D3) 
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APPENDIX E: 



We now compute the integral defined in (0). Again resorting to the change of 
variables C, u, 



which is equal to 



—36 

^2 = 7^77-^/ rfuini(l-ni)SFi(6,6 + l;6 + 2;ni) . (E2) 
¥[b + 1) JO 

This last integral can also be done and we finally obtain 

3J-^('-.'-+1.2;.+ 2,. + 3;l). (E3) 



After some rearrangement, we derive 

J2 = -36 



(E4) 



involving Riemann's zeta function ({z, q) = Z^cTl'^ + ^ ■ Plata's formula of Appendix 
may be used in the zeta function, yielding 

^ 2(6+1)2 ^ (6+1) e2-?/-l [y2 + (5 + l)2]2 ' ^^^^ 

so that the result for J2 is 

n- _ na\ ^ r°° dy 4(6+1)?/ 1 

~ 1262(6 + 1)2 "^io e2-3/-l [y2+ (5+l)2]2 j • ^ 

If we expand the last fraction in ( [E^ ) in a power series in y, use the definition of the Bernouilli 
numbers, 

B2n _ , r dyy 



ir-r^F^, (E7) 

' Jo e^^y - 1 ^ ^ 



4n 

and the fact that B2 = 1/6, we finally obtain (|53D . 

APPENDIX F: 

We compute the effective potential at zero temperature for c? = 3, in order to illustrate 
the workings of the renormalization procedure. The calculation can be viewed as a special 
case of the one performed in the text, with the kink being traded for a constant background 
ifc- In that case, determinants are simple to compute, and we formally obtain: 
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1 



A r d^k 



In 



1 + 



2(A;2 + M2) 



The subtraction of counterterms leads to: 



(Fl) 



(F2) 



which will yield the renormalized expression if we use a cutoff A to regularize (|F2|) , and fix 
renormalization conditions. Since we may write the renormalized effective potential Vr as: 



n=l 



(F3) 



it is easy to show that the renormalization conditions (|T7| ) and (|58[) determine Ci and C2: 

A /-A 1 



(27r)4 (P + M2) 



(F4) 



' 4 7 (27r)4 + M2)2 • ^ 

Clearly, both expressions correspond to one-loop graphs at zero external momenta. The 
renormalization procedure adopted in the text also uses zero-momentum subtractions in 
expansions in {^p^ — Lp1)/2, but in the kink sector. In that case, the effective action in the 
kink sector replaces the effective potential, while the expansion ( [F^ ) becomes a functional 
one 

°° 1 /■ d'^pi f d'^p ~ 



with <^(p) standing for the Fourier transform of the kink background. 
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